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J^>- Abstract 

In the article, correct method for the kinetic Boltzmann equation asymptotic solution is for- 
mulated, the Hilbert method and the Enskog error are considered. The equations system of 

<■> ' 

multi-component nonequilibrium gas-dynamics is derived, that corresponds to the first order in 

>> : 

the approximate method for solution of the kinetic Boltzmann equation within the Struminskii 



approach. It is shown, that velocity distribution functions, received by the proposed method and 
by the Enskog method within the Enskog approach to solving of kinetic Boltzmann equation for 
gas mixture, are equivalent to first infinitesimal order terms (inclusive; accordingly, systems of gas- 
£f~) , dynamic equations of the second order coincide), but, generally speaking, differ in the next order. 

o 

lO ' An interpretation of turbulent gas flow is proposed, as stratified to components gas flow, which is 

o ' 

described by the derived equations system of multi-component nonequilibrium gas-dynamics. 
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I. INTRODUCTION 

In 1912 Hilbert in [1], Chapter XXII, as an example of integral equation, considered 
the kinetic Boltzmann equation for one-component gas and proposed a "recipe" for its ap- 
proximate (asymptotic) solution. Enskog concretized the Hilbert "recipe". But meantime 
Enskog had made a not obvious logical mistake (see below section IIIip and, as the result, 
had formulated untrue method for (asymptotic) solution of the kinetic Boltzmann equation, 
having proposed to use "null" conditions, conditions (jl2l - ([Hj) below with zero right-hand 
side, for determination of five arbitrary functional parameters of the first and following ap- 
proximations of the velocity distribution function. As result of paralogism of the successive 
approximations method, partial time derivatives vanish in the necessary and sufficient con- 
ditions of solutions existence of integral equations of higher orders, see equations (H61) below, 
and with them terms of gas-dynamic equations, corresponding viscosity, heat conductivity 
vanish. Enskog "improved" the situation by the introduction (see, for example, [2j, 
Chapter 7, § 1, Section 5) of the unreasonable expansion of partial time derivative. 

Below, on the example of proposed by Struminskii in [s| method for approximate solution 
of the kinetic Boltzmann equation for multi-component gas, it is shown, how it is necessary 
to change the Enskog method. The Struminskii approach differs from the Enskog approach 
to approximate solving of the kinetic Boltzmann equation for gas mixture (see, for example, 
j^], Chapter 8, § 2) in the expansion of the Boltzmann equation - see (j2j) below. It should be 
noted, that approaches to the approximate solving of the kinetic Boltzmann equation, that 
are close to that suggested by Struminskii, were considered previously in the kinetic theory 
of plasmas, see, e.g., J], § 7.5. Struminskii's paper 3| and, hence, his following papers, 
referring to [s| (for example, Q]), involves also errors in calculation of collision integrals. 
This led Struminskii to improper conclusions of general nature. Proposed change of the 
Struminskii method does not clear up the principal demerit of the Struminskii's method, 
i.e. the lack of the explicit physical small parameter, in terms of which the expansion of the 
successive approximations method is performed (but it is not very significant for the content 
of this article, as we shall see below in the section IVj) . 

Further in this article calculations of collision integrals of the general form and for a 
specific model of the rigid sphere potential for the case, where separate components have, 
generally speaking, different mean velocities and temperatures, are presented. The equations 



2 



system of mult i- component nonequilibrium gas-dynamics is derived, that corresponds to the 
first order in the approximate method for solution of the kinetic Boltzmann equation (the 
approximation order is defined below) within the Struminskii approach. It is shown, that 
velocity distribution functions, received by the described method and by the Enskog method 
within the Enskog approach to solving of kinetic Boltzmann equation for gas mixture, are 
equivalent to first infinitesimal order terms (inclusive; accordingly, systems of gas-dynamic 
equations of the second order coincide), but, generally speaking, differ in the next order. 
This difference is the possible reason, that going to higher order in the Enskog method does 
not lead to any real improvement in the results. 

At the end of this article an interpretation of turbulent flows is proposed within the 
multi-component gas-dynamics. 

The below notation is close to the one in {2]; this allows an easy comparison of the 
described below theory to the Enskog- Chapman theory and replacement of the treatment of 
details, common to the theories, by references to the appropriate points in [2|. Also it may 
be useful a familiarity with basics of the asymptotic expansions theory, accurately stated, 
for example, in [gj], Chapter V. 



II. MODIFIED STRUMINSKII'S METHOD FOR KINETIC BOLTZMANN 
EQUATION SOLUTION 

The basic idea of the Struminskii's approach is as follows. In Boltzmann equation for 



mixture of rarefied monatomic gases [see 



Chapter 8, (1.1); the derivation of the Boltz- 



mann equation and its domain of applicability are d 1S cu SS ed, for example, in g Chapter 3 
and 18, [21], Chapter 7, § 1 and the Bogolyubov paper |8J, which is included in [2j as an 



addition! 



dt 1 dr win dc,; 



J2JJJ (nfhf^)9ijbdbded Cj 

E// (fifj-fifi)^^^, (1) 



3 



the formal parameter 9 (used simply as an "indicator of smallness", 9 = 1) is included in 
somewhat another manner, than in the Enskog approach: 

T>iU = ~\ •/,(/;./! ^- /;,(/,//!. (2) 
In the Enskog approach multiplier 1/9 refers to the whole right-hand side: 

Vifi = -- e Y,Ji3Ui,h) (3) 

n 

- see [2j, Chapter 7, § 1, Section 5. 

In ([I])-© subscripts y number components of the mixture; Xj - the external force, acting 
on molecule of i-th grade; rrii - the mass of molecule of i-th grade; gr^ - the module of relative 
velocity = Cj — Cj of colliding particles; h - aiming distance, e - azimuthal angle, k is 
unit vector directed to the center of mass of colliding particles from the point of their closest 
approach to each other, see |2j, Chapter 3, figure 3; the scalar function ky (gy,k) is defined 



by the equality 



in (H 



Qijbdbde^ kijdk.] (4) 



Ji(fij) = fj (fif-mhdkdc, (5) 

Jij (/<, fj) = J J (fifj - fifS) % dk dc 3) (6) 

to differentiate between velocities of colliding molecules of the same grade in (151) , we will 
denote one velocity by Cj and the other by c (without any subscript) and omit the subscript 
in the relevant velocity distribution function f - cf. with [2J, Chapter 3, § 5, Section 1; the 
other notation is essentially the same as in Q], Chapter 3, § 5 and Chapter 8, § 2. 

Formally introduce parameter 9 to the series of successive approximations for the velocity 
distribution function ff 

h = ff ) + e f p + 9> f p + .... (7) 



Write the differential part of equation §1§ as: 

^ r ( 9 d Xi d 

T>iU = 7^T + c i-^- + ^~ 

\ot or rrii ac, 

i AT>(2) , n2^(3) 



v w + evf 1 + 9 A vf > + ■■■ , (8) 

9t ' <9r m,- 9c,' 



0], Chapter 7, § 1, Sections 4, 5 and js]. 



- cf. with [2j, Chapter 7, § 1, Sections 4, 5 and [3j. In dHJ)-© the partial time derivative 
expansion 

(10) 

dt ^ dt y K ' 

used by Enskog and then by Struminskii, is not. As a result, the below-discussed method 
for the Boltzmann equation solution differs in essence from the Enskog's and Struminskii's 
methods. 

Having substituted ([7]) and ([8]) into ([2]) and equated coefficients in like powers of 9, we 
arrive at the system of equations of the successive approximations method for finding of 
functions f { ■ . This system of equations can be written as [cf. with [2J, Chapter 8, (2.5)- 
(2.8)]: 

Pf + E^ + ^ + ^O, (11) 



where 





= 0, 


t(0) 


= 0, 


J ij 


f(0) 


= 0, 


f(0) 


= Ji 


j(l) 


= 0, 



(12) 
(13) 
(14) 

?°¥ o) ) + ^ (/r/ (o) ) , (15) 

(16) 

4 ) = ^(^ 0) r i) )+--- + ^(^" 1) ^° ) )' ( 17 ) 
J? +1) = Ji(fPf {r) )+--- + J i (ft ) f {1) ), (18) 
J\ r) = Jih^A+Jitff'f®), (19) 



r = 1,2. . . . 

The assumption that X>j r ^ does not contain , reduces the problem of the integro- 
differential equation solution to a simpler problem of the integral equation solution. 

The formal parameter 9 is not introduced to the Boltzmann equation in the Enskoe 
approach ([2|, Chapter 7, § 1, Section 5), but the 9 is introduced to the series of successive 
approximations for the velocity distribution functions in another manner than in (J7j): 

f> = lf? ) + f? ) +0f! 2) + ---. (20) 



The result is the same, if expansion ([7]) is used for the velocity distribution function, but 
multiplier 1/8 referring to the whole right-hand side is introduced to equation ([3]). The 
successive approximations // ,// ,// ••■ calculated within the Enskog approach prove 
ordered in (inverse) mixture molecule number density n: is proportional to n, is 
independent of n, etc., - see j^j], Chapter 7, § 1, Section 5 and Chapter 7, § 2. Ipso 
facto, a physical substantiation of the using of the successive approximations method for 
finding of a solution of the kinetic Boltzmann equation appears. For Struminskii's approach 
it is inconveniently to define explicitly a small physical parameter, in terms of which the 
expansion proceeds. 

Below, when discussing the approximation order, we will consider the approximation 
order equal to the value of the superscript r in (fTTj) . In the zeroth approximation we obtain 
the following integral equation for determination of function f^: 

^(/• (0) ,/ (0) )=0. (21) 

Having multiplied equation ( l2Ti) by In /», integrated over all Cj, and transformed the integrals 
with taking into account 





flf'k ij dkdc i dc j = /// (j)' i f i f j k ij d\iidc i dcj, (22) 



with equality (122]) being independent of the form of functions <f), f (cf. with [2], Chapter 3 



3, Section 3), similarly to [2], Chapter 4, § 3, we find: 

j(0) yr(O) 




f jf» f W_ f jfS) f V) ) kidkdcdc . = . (23) 



y(0)'y(o)/ 

The integrand in (1231) cannot be greater than zero, therefore the integral in (1231) can be equal 
to zero only when the integrand (it is supposed, that all considered functions are continuous 
in each point of the range of their definition) vanishes for all values of the variables, i.e. 

/f 7 (0) ' = /f 7 (0) (24) 



or 



In + In /<°> - In - In /<°>' = 0. (25) 
Hence, In fj; should be expressed linearly in terms of additive collision invariants ip^ = rrii, 

I (2) ; (3) 1 



2 

2 ""i^i ■ 
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For the collision of the i-th molecule with the j-th molecule the conservation of invariant 
ipf^ is expressed by equality: 

4 l >' + V>f - 4 l ) - = (* = 1,2,3). (26) 

As molecule velocities c^, c'j upon the collision are determined completely in terms of the 
molecule velocities c i; Cj before the collision by two parameters specifying the collision, for 
example, by aiming distance b and azimuthal angle e (see above), from four scalar equations 
[261) . corresponding to conservation of energy and three components of momentum, like in 
2J, Chapter 3, § 2 we conclude, that no additive collision invariants linearly independent of 
ipf exist. However, this conclusion seems not quite stringent logically. 
Thus, 

ln/f = af< 0) + ■ m lCl - af^-m^l (27) 

where and a- 3,0 " 1 are some scalar functions of r and t, independent of c i} and ct[ 2 ' ^ is 

a vector function of r and t. Or 

where is a new scalar function of r and t. I.e. the general solution of equation ( |2T|) 

can be written in the form of Maxwell function: 

fr = pr H^\ i e -^''YH'\ (29) 



2nkf3, 



(3) 

i 



where 



A W = «. M HpT ■ (30) 

(2,0) 

A (2) = %sr, (31) 



( . ( :,m 

}(3) 



By definition particle number density, mean velocity and temperature of the z-th compo- 
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nent are introduced: 

ni = f J fidci, (33) 

n, m ,u, S /m.c/.dc,, (34) 

^ n^Tj d = J ~ m, - Ui) 2 £ dq. (35) 

From (|33l -(l35l) we obtain equality: 

^ n^T- + i n^w 2 = J ^ miC 2 /j dcj, (36) 

which is convenient for using further instead of definition fl35l) . 

Together with asymptotic expansion ([7]), according to definitions (|33|) . (|34|) . (|36|) . it is 
necessary to define asymptotic expansion for particle number density rii of the i-th compo- 
nent 

n . = n f) + 0n W + 2 n (a) + . . . ; (37) 

mean velocity Uj of the i-th component 

u, = uf + OvP + 9 2 uf } + • • • (38) 
and temperature Tj of the i-th component 

Ti = rf } + er^ + e 2 Tl 2) + ■■■ . (39) 

Having substituted ([7j) and (13"T|) - (l3"9"j) into (1331 . ([3"I]) . (13"B1) and equated zeroth infinitesimal 
order terms, arbitrary functions fi^ (r, t), (3^ (r, t) and (3^ (r,t), that appear in (129|) . can 
be juxtaposed with the zeroth approximations to local values of particle number density, 
mean velocity, and temperature of the i-th component: 

A (1) M) = nf } (r,t), (40a) 
/3S 2) (r,t) = uS 0) (r,t), (40b) 
A (3) M) = T/ 0) M). (40c) 

Functions (l4"0"j) are found from the first-order gas-dynamic equations system (below we 
will see, why from this equations system): 

/ ^ (V } + E J S ] ^j dc * = ° ( z = !> 2 ' 3 )- ( 41 ) 



Similarly, functions (r = 1,2. . .), found at the r-th step of the successive approxi- 
mations method, prove, see below, parametrically dependent on 5 arbitrary scalar functions 
of r, t. Having equated terms of the same infinitesimal order in (1331) . (jMj) . (1361) : 



frdci = nr, (42) 
/ rriiCifP dci = rrii (njUj) (r) = m, E 77 ^ ~ s) uf ) , (43) 

^ s=0 



= 2*E + i m, £ E n<-'>u<-*> • u®, (44) 

s=0 ' s=0 q=0 

we obtain 5 conditions for each subscript i, that can be used to express 5 arbitrary functional 
parameters in / 4 as functions of 

(45a) 
(45b) 
(45c) 

Functions ( 1451) can be found from the (r + l)-th order equations system: 

/ ^ + E + 4 r+1) ) = (* = 1. 2 > 3)- (46) 

In ( l4Tj) and (J46]), to simplify the further transformations according to the definitions of 
the pressure tensor of the i-th component 

Pi = f J rrii (ci - Uj) - Ui) /< dc, t 

= J rriiCiCifi dCi - r^m^U;, (47) 
the vector of heat flow density of the i-th component 

<«-«.)' <«-,.,)/,*, 

/l 3 
- rriicfcifi - p< • Ui - - riikTiUi (48) 

and the temperature of the i-th component (|35|) . \E , J - 1 ' 1 = m*, ^ = m^Ci, = ^rriiCf, 
where Cj = (c^ — u$), can be used instead of ipf\ 



(r) 


= n^(r,t), 


(r) 
u i 


= u^M), 


T (r) 
i 


= T«(r,t), 
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Functions f$ , r = 1, 2 . . . are solutions to integral equations (fTTj) . which can be rewritten 

as 

Vr ) + Y,4 ) + Jt ) = -Jt ) - (49) 

The left-hand side of (j49p includes only functions known from the previous step of the 
successive approximations method. Unknown function f> r ^ appears linearly only in the right- 
hand side of equation (j4"5|) . Therefore, the general solution of equation fH9|) is + Q r \ 
where is the general solution of homogeneous integral equation 

4 r) = ^(/f/ (r) )+^(/r/ (0) )=0, (50) 

and is some partial solution of inhomogeneous equation (|49|) . 

Solutions Hj and £^ will be sought in the form of E^ = f^^>^ and £^ = fj; ^(j)\ r \ 
where ^ and 0- are new unknown functions. In view of (124|) . 

= // if V (0) + 0S r) - W ' - <AP) feidkdc 

d = f n?I 4 = 0. (51) 

Having multiplied equation ( 1511) by 4 dc,, integrated over all Cj, and transformed the 
integrals with account for (|22|) . we arrive at: 

\fjf if V (0) (> + # - - 2 rfk dc dc, = 0. (52) 

From (1521) we conclude, cf. with (123]) and (127|) . that <p( is a linear combination of additive 
collision invariants ipf . 

At) (l.r) ; (1) . (2,r) ,(2) . (3,r) / (3) /r- \ 

where cq and of are arbitrary scalar functions of r and t, and of is an arbitrary 
vector function of r and t. In place of the additive invariants \ functions can also be 
used, which are additive invariants of collision of particles of the same grade: 

4 ) = ^ r) '^+^.^+^?\ (54) 
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where aj 1 '^' and af^' are new arbitrary scalar functions of r and t, and ot^ 2 '^' is a new 
arbitrary vector function of r and t. Thus, 

^ = f^(4' r) + ^ r) -mc i + at ) lrn^,). (55) 



To make use the results of the integral equation theory, transform equation (j49[) to the 
standard form. The right-hand side of integral equation (149j) . i.e. nfli (^^ r '- ) ), which is a 
function of Cj (and, naturally, of r and t, to simplify the notation, the evident dependencies 
are not specified) can be represented as 

n?Ji(#M) = K (a) ^ (Ci) 

+ y if (c^, c) <Z> (r) (c) dc, (56) 



JJfPf^dkdc 

fP [[[ f (0) \^-c\bdbdedc, (57) 



where 

K ( Ci ) 



and if (cj, c) is a symmetric function of Cj, c (see Chapter XXII and [2], Chapter 7, § 6). 
Hence, equation (jl9|) can be rewritten as: 

if) = if (c,) <*f (c) 

+ y if (ci, c) $W (c) dc, (58) 

if (cj) in ( J58l) denotes the left-hand side of integral equation ()49l) . Linear integral equation 
( 1581) is reduced by transformation 

t \ $ {r) (c-) 

*, W (c*) = 4=, (59a) 
V (Ci) 

if (c) = if (c,) v^oM, (59b) 
if( Ci ,c) = if (ci,c) V^o(ci)if (c), (59c) 

maintaining the kernel symmetry to the linear integral equation of the second kind with 
symmetric kernel 

if } (c^) = W (a) + y if (a, c) *w (c) dc. (60) 
n 



As homogeneous integral equation (15T|) has nonzero solutions ( |53|) . corresponding to equation 
fl60|) . homogeneous integral equation 



= U r) fc 



+ y ^( Ci ,c)0M(c)rfc. (61) 



has nonzero solutions 



x 



V^M^)- (62) 



Therefore, according to the second Fredholm alternative ([9J, Chapter III, §§ 2, 3 or 10], 
Section 15.3-7), if K (cj, c) is piecewise-continuous and normalizable and (cj) is contin- 
uous and square-integrable (these conditions are assumed met), the necessary and sufficient 
condition of solution existence of equation (160]) . with taking into account symmetry of kernel 
K(ci,c), is the orthogonality of (cj) (with weight function 1) to each solution 0^(cj) 
of equation (IBTjl . 

Thus, the necessary and sufficient condition of the existence of a solution to equation 
is the orthogonality of the left-hand side of the equation to functions ip^ (c$) [or, which is 
the same, to functions ( c i)]j i- e - satisfaction of equalities (jUJ, (1161 . which, on the other 



hand, can be considered as equations for determination of n\ , uf \ , r = 0, 1, 2 . . .. 

The partial solution of inhomogeneous integral equation flU sf } = // 0) ^ r) can be 
constructed, for example, using expansion #j(r) (cj) in series in terms of Sonin polynomials 



with expansion coefficients, depending on r and £, as this is done in {2], Chapter 7 and 8 



III. SOME REMARKS 



By the successive approximations method we receive, generally speaking, asymptotic so- 
lution of the task. Upper statement of the method for the kinetic Boltzmann equation 
solution is close to [2|, Chapter 7 and 8. Somewhat more accurately the method for the 
kinetic Boltzmann equation solution can be described within the theory of asymptotic ex- 
pansions with variable coefficients [6], Chapter V, § 2, Section 5. Additional conditions (for 
example, initial or boundary conditions for a differential equation), if there are ones, it is 
necessary to expand into asymptotic series, and to solve, equating coefficients at same terms 
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of unified scale of comparison - Chapter V, § 2, Section 1 (usually, the scale of compar- 
ison is the set of functions 8 r , defined on the neighborhood filter of the point 6 = 0), the 
resulting system of equations. If this system of equations can be resolved, we have asymp- 
totic solution of the starting equation, satisfying supplementary conditions. Sometimes the 
found asymptotic solution is regular Chapter 1), i.e. analytically depending on 9, 



solution of the task. For example, the equation (see [12| , Chapter V, § 2 or [13] , Chapter IV, 
§7.1) 

of f +f = 0, /(0) = (63) 

with zero initial condition have with regard to comparison scale of function 9 r , asymptotic 
solution (t) =0, (r = 0, 1, 2 . . .), that is also the exact solution of the task. However, 
the task 

of f +f = 0, /(0) = 1 (64) 

has not analogous asymptotic solution, as the zero-order asymptotic solution of the differen- 
tial equation (t) = contradicts to the initial condition (0) = 1. This is not serious 
weakness of the successive approximations method. In the task f[M|) one can introduce new 
function g (t) : 

f(t) = exp{-t/0)g(t) (65) 



cf. with ll|; though discussed by Lomov expansions are not asymptotic expansions with 



variable coefficients, and accordingly Lomov' s approach as a whole seems not quite correct, in 



the monography [llj is actually shown, that similar replacements of functions allow to obtain 
asymptotic solutions for a wide class of problems) and by the successive approximations 
method receive asymptotic solution of the task g^ (t) = const = 1, g( r > (t) = const = 
(r = 1, 2 . . .), that is again the exact solution. 

The criticism of the successive approximations method [12j, Chapter V, § 2 or [13), Chap- 
ter IV, § 7.1, possibly, reflects a dissatisfaction of authors with the unreasonable expansion 
of partial time derivative in the Enskog method. 

Hilbert, having marked in [1], Chapter XXII, that the expansion 

F = ^ + $ + XA + ---, (66) 
A 
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analogous to (1201) upper, (Hilbert considered only one-component gas; we maintain here 
Hilbert notation, however clear enough from a context) is the power series in (small param- 
eter) A, satisfying to the Boltzmann equation and such, that expressions [cf. with (13"31- (13"6"|) 

and fl37D-fl3HD] 

+ A J \p (i) Xduj + --- (i = 1,2,3,4,5) (67) 

for t = to pass into power series 

f(0 

A «) = J +g ® + X h® + -.. (i = 1,2,3,4,5), (68) 
A 

in the theorem, closing his work, has formulated "recipe" for (asymptotic) solution of the 
kinetic Boltzmann equation, in which he has proposed five arbitrary functional parameters 
of functions . . to define "from five partial differential equations" , analogous (|4ip . 

(HEP , " at that for t = t " to preset 

J i) {i) $duj = AA W (i = 1,2,3,4,5), (69) 
J^ i) ^duj = (i = 1,2,3,4,5), (70) 

iP {i) Xduj = (i = 1,2,3,4,5). (71) 



In notation from fl33|) - fl36|) . (j3"T]) - (j3"9~]) . (JH1) . (H6|) upper, Hilbert proposed simply to use special 

initial condition 

n(r,t ) = n^(r,to), (72) 
u(r,t ) = u(°)(r,t ), (73) 
T(r,t ) = T(°)(r,t ) (74) 



or 



= (/= 1,2,3) (75) 

t=to 



for r = 1, 2 . . .. It allowed him, as corollary of the theorem, to formulate "fundamental result 
for the theory of gases: the state of stable gas at any t is uniquely determinated, if for it at 
t = t density, temperature and velocity are known as function of a point of space" . 

14 



"For the further substantiation of the gas theory" it would be necessary to supplement 
Hilbert's theorem with explicit definition of five arbitrary functional parameters of functions 
f\ r \ found on the r-th step (r = 0, 1, 2 . . .) of the successive approximations method, through 
gas physical parameters f[3"3"j) . f[3"4"l) . (j3"5]) . (j4"2"|) - flSJ), but Hilbert had not made it. 

Enskog formulated up the Hilbert "recipe" for concrete calculations. However meantime 
Enskog had made a logical mistake. He used "null" conditions ( IT5]) identically, at any t, not 
just at t = to (see Chapter 7, § 1, Section 1): 

J^f^dc = (1 = 1,2,3) (76) 

for r = 1,2.... From the point of the successive approximations method view Enskog instead 
of fl37|) - fl39|) had supposed 

n(r,t,9) = 6° n (r, t, 6) + 6 1 + 6 2 + • • • , (77) 
u (r, t, B) = 6° u (r, t, 6) + 1 + 6 2 + ■ • • , (78) 
T (r, t, 9) = 6° T (r, t, 6) + 1 + fl 2 + ■ • • . (79) 

If n, u and T did not depend from r and t, it would mean, that Enskog used simultaneously 
different scales of comparison {n (6) , 6\ 6 2 . . .}, {u (6) ,6\6 2 . . .}, {T (6) , 6\ 6 2 . . .} in the 
successive approximations method, that is already wrong. In a general case, when n, u 
and T depend from r and t, the sums (J7"7"]) - (JT9"]) cannot even be considered as asymptotic 
expansions with variable coefficients. 

Infringement of logic of the successive approximations method is immediately appeared 
in that from the equations, analogous (1461) (r = 1, 2, . . .), in compliance with (|76ll partial 
time derivatives vanish 

J ^ {l)dJ ^dc = ^J ^f^dc = (/ = 1,2,3), (80) 

and with them terms of gas-dynamic equations, corresponding viscosity, heat conductivity 
. . . , vanish. Somehow to correct the situation, Enskog has been forced to enter unreasonable 
expansion of partial time derivative ffTUj) . 

IV. CALCULATION OF DEFINITE MULTIDIMENSIONAL INTEGRALS 

In this section we are dealing with calculation of definite multidimensional integrals 

//// ^ ( /;(0) ^ 0) " tobdbdedad*. (81) 
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In (EH) *W = m i? *f> = m^Q, ^ (3) = fr^C?, d = d - u, 



3 



// 0) = ^(^) 3 e- mi(Cl - Ul)2 / 2fcr % (82) 

is the Maxwell function of distribution of velocities of the z-th component particles, the 
prime in the distribution function means, that the distribution of the particle velocities d i 
after the collision is considered. The other notation is specified above. 
According to (122]) . integral flHTj) can be transformed as follows: 




*i° (f^f^-f^f^g^dbdedc^ 
^f^f^g^b'db'd^d^ 



^ i i ) ff ) ff ) g t] bdbdedc i d Cj 






- /f db de d Cl d Cj . (83) 

As the particle mass is conserved in the collision, for = rrii integral (1831) vanishes. 
In the two other instances, generally speaking, this is not the case because there is no 
summation over the components, cf. with p|, Chapter 7, (2.33). 

Hereafter statements of the two following simple propositions are used several times. 

Proposition 1. / is assumed to be a ruled function on R with values in R, w £ R 3 be a 

fixed nonzero vector, n £ R 3 be a unit vector. In this case 

7T 

J f (w n) n dQ a = ^ J f{w cos (0)) cos (0) sin (9) dO. (84) 



o 



In the left-hand side of |jO ) the integral is taken over all directions of vector n, w • n is the 
scalar product of vectors w and n. 

Remark. If w is a zero vector, then the right-hand side of (8$ ) is set equal to 0. 

Proof. Select the system of spherical coordinates, such that the polar axis direction be the 
same as the direction of the vector w. Resolve the vector n into two components: parallel 
(ny) and perpendicular (n_|_) to the vector w: 

(w ■ n) w . 
n = ni, + n ± = ^ J— + n ± . 85 
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Having substituted expression (1831) for the vector n into the left-hand side of (184)) and inte- 
grated over the azimuthal angle, we obtain the required equality ({83]) . as in the integration 
over the azimuthal angle the containing term vanishes. □ 

Proposition 2. E and F is assumed to be two complete normalized spaces over field R 7 u 
be a continuous linear map of E into F. In this case, if f is a ruled function on interval 
I C R with its values in E, then u o / is the ruled function on I with its values in F and 

Ju(f (t))dt = u [Ji (t)dt J . (86) 



Proof. Equality (1861) follows immediately from the expression for the derivative of composite 
function uo/; the details of the proof can be found in [6j, Chapter II, § 1, Section 5. □ 

In these propositions ruled functions can be replaced by better known continuous func- 
tions. 

The major difficulties in the calculation of integral ( 1831) are associated with the fact that 
parameters of the Maxwell functions for the i-th and the j-th components are not equal: 

\ii ^ uj, Ti ^ Tj. (87) 

As a result, it is not easy get rid of the scalar products of vectors in the exponent (it is 
desirable that the expression for the exponent be as simple as possible). 

As the scattering angle depends on the module of relative velocity of colliding particles 
[see, for example, [2], Chapter 3, § 4, Section 2 or Q], Chapter 1, (5.26)], it is natural to 
transfer in (1831) to new variables - center-of-mass velocity Gy and relative colliding particle 
velocity gy, which are related with the particle velocities and Cj as: 

d = Gij + ^ — 88 

rrii + rrij 

Tfl ■ 

c 3 = G ij T — Sij, (89) 

VTLi + TTlj 



a, 



- cf. with [2J, Chapter 9, § 2. For further simplification of the exponent vector Gy can be 
replaced by vector G^ resulting from G^ in an arbitrary affine transformation, for example, 
the one, which is a composition of shift, homothety (multiplication by a scalar), and rotation. 
The rotation arbitrariness is reduced to the freedom in choosing of direction of the polar 
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axis in the transition to the spherical coordinate system. Similarly, the vector can be 
replaced by the vector g^, resulting from g^ in composition of arbitrary homothety and 
arbitrary rotation. The shift of the origin of the vector g^ would lead to a parametric 
dependence of the final integral on vectors Uj and uj (cf. with [14j . Chapter 3), which is 
undesirable, as integral (1551) is supposed to be reduced to Chapman-Cowling integral fi£ 
[see fl, Chapter 9, § 3, (3.29) and Q, Chapter 7, (4.34)]. 

In view of the aforesaid, make the following substitution of variables G^- and g^-: 



5'J 



z\ Bij, (90) 
G tj = z 2 Gij + z 3 g tj + Ul + Uj . (91) 

In ( 1901) - (19T1) the scalar factors zi, z 2 , and Z3 are selected from the condition that the coeffi- 
cients of gfj and Gjj in the exponent be equal to 1 and the coefficient of the scalar product 
gij ■ be equal to (compare to the method of variable separation): 



z x = iHmTj + mjTi)^ (92) 



VTLiVTLj 



Z2 = J T , (93) 

2s = ^r|^z. (94) 

rrii + rrij y 2 [rriilj + rrijli) 

Analogous substitutions of variables can be used in more complicated situations, for example, 
discussed in [14j . Chapter 3. 

With new variables the exponent can be written in the following form: 



+ G\a + a w 2 + aigjj • w + a 2 G,j • w , (95) 



where 



w 



Uj — Uj 

2 ' 

; rr, 

2T] + 2Tj 



(96) 



rrii rriA 

a = + -f , (97) 



2m,mj 

ax = -2 X ^ 98 

w rriilj + rrijli 



, rrij rrii\ / 2T{Tj 
Tj Ti J V rriiTj + rrijTi' 1 
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It is easy to see that it will be impossible to get rid of the constant term in exponent (|95|) 
and, hence, of the constant exponential factor, which will appear hereafter in all expressions 
containing integrals of form (1ST]) . (I8"3"j) using only the above-specified transformations of 
variables (without using the shift of the origin of vector gy). Such factors are missing in jjj, 
(8). 

Determine Jacobian of transformation of variables (cy cy) — > (gy,Gy ] [see ( )90]) - (l9Tj) ]: 



d(Ci,Cj) d(Ci,Cj) <9(gy, Gy) _ 3 3 d(Ci,Cj) 

Z l Z 2 



d (gy, Gy) 8 (&> G ^ 8 (gy, Gy) 9 (gy, Cj + ^ gy 

_ r 3 „3 d (Cj, Cj) _ 3 3 d(Cj,Cj) _ 3 r 3 M nn\ 

~ Z ^d (gy, Cj ) " Z > Z * d (C t - Cj , C,) " ^ (1UU) 

Now consider the case, where *f = *f } = m, (q - u^). In view of ([95]), (jlODjl . (|90|)-(l9T|) 
and the equality, following from the definition of k upper, 



rm « - *) = (gy - gy) = -2 -^p- ( g .. . k ) k (101) 



integral (1831) can be rewritten as: 




mi (c- - a) fi 0) fj 0) 9ijb db de dc { dcj 



-2 mim i ztzlnA^-) nA^-\ 1 1 1 (gy - k) k 
rm + rrij 1 2 \2nkTiJ 3 \2nkTj 1 111 VS j ; 



3 3 
2 / <m ■ \ 2 




x exp I — 



9ii + + a w 2 + aigy • w + a 2 Gy ■ w 



x ^y 6 de e£Gy dgy <i&. (102) 

Integrating with respect to e in (11021) (with fixed gy and Gy), resolve vector k into two 
components: the ones parallel and perpendicular to vector gy - cf. with the proof of 
Proposition CD 

J (gy ■ k) kde = 2ttcos 2 - gy = 7T (1 - cosx) gy- (103) 

When integrating over Gy and directions of vector gy, use Proposition [TJ As a result we 
arrive at 

uii (c ■ - Cj) // 0) fj 0) gijb db de dc { dcj 




rriiTj + rrijTi w ^/ti 



,2 



16n,-n,- V — e m * T i+ m J T » 



* J (m ; + 771,-) w £ 2 



x j j e 9 *j [giji cosh (#y£) - sinh (#y£)] #y (1 - cos x) b db dg ip (104) 
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In (TT04D 



(105) 



factor a\ is determined by formula (|98|) . It is easy to check that the singularity at £ = 0, 
which is possible when if = 0, is actually absent in the right-hand side. Expression fjl04j) 
differs from Sruminskii's expression, i.e. |3|, (8). 



The case, where ^ 



(0 



(3) 



(cj — Uj) 2 , differs from the just considered one in 
the factor of the exponent in the right-hand side of (11021) . Transform difference — 
according to (1901) . (19T|) and Chapter 3, (4.9) and taking into account that only the 
relative particle velocity direction changes during the collision (g^ = g'^): 



(3)/ 



^(3) = mt 



vrtiVrij 
rrii + rrij 



U; 



- C; 



U; 



(c- + Cj - 2Ui) 

({^. - gy } • {Gtf - U,}) 



"2 2l (gy • k) 

rrii + m,- 



X ( k • { Z 2 Gij ; + Z 3 gij 



11 - u, 



(106) 



With respect to its arguments the scalar product is a bilinear continuous function, there- 
fore Proposition [2] can be applied. On the integration with respect to e, similarly to (11031) . 
we arrive at: 



-2z 1 



mi + m 



k) k • \z 2 + z 3 - 



U; — U, 



rriimj 

-2n Zi (1 — cos x) 



%2 Gij + z 3 gij — 



de 

Uj — Uj 



rrii + rrij v "" V~ J I " " J "~" J 2 

Perform the integration over G^ and directions of vector gij using Proposition [T] 




rrij 
2 



Ci - u, 



fi^ff^9ijb db de dci dcj 



16nj?2j^— e 



-5 2 . 

e y ^ 



w 



x < - cosh (<7y£) - sinh (<fo£)] - ID^g^ sinh (<fo£) 

^2 



Xfe (1 - cos x)bdbdgij. 



(107) 



(108) 
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In (TT08D : 



Di, 



2 rrijTi 



(109) 



rrii + mi ' 



D 



miirtj (Ti — Tj) 2 Ti 2 Tj 



(110) 



2 (mi + m 3 f V m rrij 



The other notations are the same as in (I104p . 

It is interesting to note that for Uj = Uj integral H104I) and the first term in (1108!) vanish 
and the second term in (11081) is: 



that corresponds to energy transfer from the "hot" components to the "cold", see the gas- 
dynamic equations system below. In view of the sign of ai ( 198]) and definition of £ (11051) . 
the first term leads to temperature increase with w^0. 

V. FIRST-ORDER EQUATIONS SYSTEM OF MULTI-COMPONENT 
NONEQUILIBRIUM GAS-DYNAMICS 

Above the gas-dynamic equations system has been derived, in a sense, as a "by-product" 
during the Boltzmann equation solution. More generally, the gas-dvnarnic eolations system 
can be written in the fan of the transport ecnfaons, of. with ,2], Chapter 3, (1.12) and 
[7j, Chapter 7, (2.31). 

Having multiplied the Boltzmann equation for the z-th component (Tj[|) by tyf' and inte- 
grated over all values of Cj (it is assumed that all the integrals obtained below converge and 
products like ^/f^Xj/j tend to zero, when Cj tends to infinity), we arrive at: 



const (Tj — Tj) , 



(111) 




(112) 
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The terms in the left-hand side of this equation can be transformed: 



dt dt 1 J dt 



( nM ] ) g^(0 



/ 



dt 1 dt 



n*-5T- > ( 113 ) 



d rW/ 

— ■n^f ) c l -n t c l -— L- , (114) 
or or 



.M^.__ffi. lT ,« X ' 



1 m, ; dd % J \ <k; 1 m, " ' 



In (1113j) - (jll5p . the bar, as usually, denotes the average of the quantity 



Vi=- I V t f t dc i; (116) 



11;, 



r and Z\ are considered as independent variables. In view of (I113p - (lll5p . from (11121) we 
obtain 



d[n&f\ « —r— /WO WO 



9 7777- dm l) d t(i) X, 



= E/// WW - gabdbdedCidcj, (117) 

- the transfer equation for the ^ , that refers to particles of the i-th grade. 

To derive the equations of mass, momentum and energy transport for the i-th component 
from flll7p . sequentially substitute m«, m» (cj — Uj), | m» (c, — u.j) 2 for ^i(Z) into (11171) . 

In the Enskog-Chapman theory, in view of the additional summation over i, the right- 
hand side of (11171) vanishes always. However, if the velocity distribution functions for some 
components are Maxwell functions fl82l) with different parameters of mean velocity and 
temperature (u.j ^ Uj, Tj ^ Tj), for example, due to some external effects (see below), 
then nonzero terms remain in the right-hand side of (I117p . In this case equations (1117)) 
(I = 1,2,3) are the same as equations (I4ip . Thus, on straightforward transformations we 
arrive at the following gas-dynamic equations system [cf. with 3|, Chapter 7, (2.42), (2.45), 
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(2.47)]: 













■ riiiii, 










~dt 


dr 






dui d 
dt dr 


(0) 

•Pi 


-Vi (0) 


= UiXi 


- riifriiUi ■ 


d 

OT 


dEi 
dt 


d (o) . (o) 


dui 
dr 


/ j e,y 
3& 


d 
dr 


■ EiMi. 





(118) 
(119) 

(120) 



In fim- (fl20]) : 



p. (0) = mmdci - u { ) (d - Uj) (0) = riikTiV = p/ 0) U (121) 



is the i-ih component pressure tensor, pf^ is the hydrostatic pressure, U is the unit tensor, 
double product of two second rank tensors w and w' ([2], Chapter 1, § 3) is the scalar 



W : W ' = Ea E/3 W^W^ = W' : W, 



q/ 0) = \ - u,;) 2 (c; - Ui) ( ) = (122) 

is the i-th component heat flux density vector, 

1 2<W 

Ei = - riim^Ci - (123) 

is the internal energy of particles of the z'-th component per unit volume, which is equal, 
in this case, to energy of their translational motion, however, the energy transfer equation, 
written in form (11201) . apparently, can be used in more general cases as well (cf. with 
Chapter 7, § 6); in fll21l) - fll23p superscript ^ denotes averaging (jllfip with Maxwell function 
if } from flE2D. 

In frrT9l-(fT2T?D 1^., J e (0 A denote integrals ffTUIj) and flTUgj) . respectively. When averaging 
the last term in the left-hand side of (I117p . external force Xj, acting on the particle of the 
z-th grade, is assumed independent of the particle velocity. 



VI. VALUES OF KINETIC INTEGRALS FOR INTERACTION POTENTIAL OF 
RIGID SPHERES 

Integral terms if^-, appearing in multi-component gas-dynamics equations system 
(I118I) -( |T20|) . are quite complex functions of mean velocities and temperatures of separate 
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components, mainly, because of a complex dependence of deflection angle x on relative 
velocity of colliding particles [cf. with [7], Chapter 1, (5.26)]: 



x(b l9 )=TT-2b [ - 1 == 



dr/r 2 



p(r) _ b 2 

w 71 



(124) 



The component temperatures appear in the resultant expressions as a making non- 
dimensional factor. 

In the simplest case of particles, interacting according to the law of rigid spheres, the 
following analytical expressions for Ip°jL, and J e have been derived: 

(o) rnjTj + rrijTi w 2 

L P ,ij ~ n i n i ,„ oc2 a ij 



X 



rrii + rrij w 2£ 



exp - 



r(0) 



2£ (e + 2) + v^F (e 4 + 4£ 2 - 4) erf (0 
4 exp (-^ [2D ltt M (e + 2) + 2D 2>ij £ 2 (£ 2 + 10)] 



(125) 



7T 



x (£ 4 + 4£ 2 - 4) + D 2iij £ (£ 4 + 12£ 2 + 12)] erf (^j 

In (TI25D-(II26D notations from ([96]), ([98]), ffT05]) . ffT09l- (fTT0]) are used. 



(126) 



VII. PROPOSED METHOD AND ENSKOG-CHAPMAN THEORY 

Let's consider, what do nonzero conditions fl42l) - fj44l) and expansions (I3~7"j) - (l39~l) for one- 
component gas lead to; for one-component gas the Enskog approach and the Struminskii 
approach coincide. 

Differential equations (141|) . from which functions n^°'(r, t), u^(r, t), T^(r,t) are found, 
for one-component gas can be written in the form [cf. with flll8p - fll23l) ]: 



n (») m ^ = „(»)X-^, (128) 
Dt or 1 

= —fE. . u (0) = _±L . u (o) ( 12 Q) 

Dt 3n^kdr ~ 3 dr { ' 
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In (H27D-(H29D 



P 



Dt dt dr' 

(o) = n (o) kT (o)_ 



(130) 
(131) 



From (I127p . fl 1 2 9 [) we receive, that in the first order gas-dynamic flow is adiabatic: 



D 
~D~t 



n 



(0) 



( T (o))" 



(132) 



For one-component gas integral equation (|49|) . r = 1, from which = f(°)$( 1 ) is 
found, with taking into account (129]) . fl56l) and d 1 2 71) - (1 1 32 j) . can be written as (cf. with 
Chapter 7, § 3): 



rr 



, J(!P a) ) = 3^ + c ^ (0) ■ x ^ (0) 



<9f 

/( o) 



<9r m 9c 

C ~ o + 2CC: — u(°' 

2 / ar ar 



In (H33D C = c - u 



(o) 



C is the module of vector C; for arbitrary second-rank tensor w 

1 . . 

w = w — — U (U : w) 



(133) 



(134) 



(135) 



- tensor with zero trace. 

The general (scalar) solution of equation (11331) . being the sum of some partial solution of 
equation (11331) and general solution of homogeneous equation / (0^) = 0, cf. with (|54"|) . we 
can seek in the form of: 



/ 



(i) 



/ 



(0) 



77,(0) y m J dr n(°) dr 



+ /(°) (V 1 ' 1 ) + ■ mC + a^^mcj 



where the vector function A is partial solution of the equation 

n/(A) = /(°) C, 
and the tensor function B is partial solution of the equation 

n/ (B) = 2/WCC. 



(136) 



(137) 



(138) 
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Solvability conditions of equations (11371) and (11381) are satisfied, i.e. (see the section HT1 
upper): 



J ^/ (O) (c 2 -0Crfc = (/ = 1,2,3), 
¥ l) f i0) CCdc = (/ = 1,2,3). 



(139) 
(140) 

In (DSHD and flUD} = m, * (2) = mC, ^ = \ mC 2 . 

Because r, t and u^°) do not explicitly occur in equations (I137f) and (11381) . and right-hand 
side of equation (11381) is symmetric tensor with zero trace, solutions A and B can be sought 
in the form of: 

A = A(n i0 \C,T i0 ~>)C, (141) 
B = B (n (0) ,C,T (0) )C°C, (142) 

where A (pS Q \ C, T"(°)j and B (n^°\ C, T^) - scalar functions of n^°\ C and T^°\ It is possible 
to impose an additional condition on the solution A (cf. with [2[, Chapter 7, § 3, Section 1): 

C 2 / (0) A(nW C,T<°>) rfc = . ( 143 ) 



Having substituted (TT36I) in (jl2l-(pi) with taking into account (j25jl. (1101 and (HHl . we 
arrive at: 



n 



(i) 



n (0) a (l,l) +mn (0) u (0). a (2,l) 



H — n 



(0) 



(3,1) 



mmu 



3A;T(°)+m( M (°)) 

mc/ (0) $ (1) dc 
mn (0) u (0) (a (y) + mu(°' ■ a^) + mAT^W 2 ' 1 ) 



(144) 



(3,1) 



(145) 



- k (nT) ' + — m(m/ 2 ) 



2\(1) 



2 J 



2 U 
1 



(0) 



3fcT(°>+m(^°>) : 



ft 



(i,i) 



H — mn' ' 

+ 5(o) kT (0) \ 3kT (o) + 2m (u^f 
4 L v ' . 



a 



(3,1) 



1 



(°)«A,(0)4 (3,1) 



+ - n w m"w w ft 
4 



(146) 
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In (j!44p - (jl46p vanishing integrals (cf. with [2], Chapter 7, § 3, Section 1) are neglected. 
From ([TllD-flUnD we have: 

n (D 3T (D 



a 



(2,1) 



n(°) 2 T(°) ' 



a 



(3,1) 



(147) 
(148) 
(149) 



n 



To first infinitesimal order terms (see. [6j, Chapter V, § 2, definition 2) expression 



/ (0) + / 



(o) 



n 



(i) 3^(i) 



a) i r wi 



u 



(150) 



n (o) 2T(°) fcT(°) A;T(°)T(°)2 
coincides with the asymptotic expansion of the solution in the Enskog- Chapman theory 



/ 



(0) 



n 



[i] 



2 — mic— u 

e V 



[ip 



'2fcTl 



(151) 



27rfcTW . 

where nM = n^ ) + nW, = u(°) + uW and = + T^\ cf. with Taylor expansion of 
function /(°) about the point (n^°\ vS°\T^°n. This assertion can be written in the form of: 



/(0) * fW + / 



(o) 



f(0) 



(i) 3^(1) 



+ t-z-ttt • mC + , „,„ N - mC 2 



u 



(152) 



y°) 2T(°) jfeT(°) fcT(°)T(°)2 

Equations (I137j) and (11381) differ from analogous equations [2|, Chapter 7, (3.9) and (3.10) 
in the Enskog- Chapman theory only in use nS \ and instead of n, u and T (i.e. n' 1 ', 
u' 1 ' and T^). In first and second terms in the right-hand side of ( 11361) . in ( 114ip and ( 1142[) 
functions nW, u(°) and T (0 \ as upper in (jSj), can be, respectively, replaced by functions 
u' 1 ' and T^. Therefore to first infinitesimal order terms expression 



/ 



(o) 



1 [2kT^\ 91nT(°) 
- A- 



77,(0) y m J dr 
coincides with the solution p 1 ' in the Enskog- Chapman theory 



(153) 



= / 



(0) 



n 



1 (2kTN\ ~ ainTW 



1 



<9 



,[i] 



(154) 



dr m 1 ' dr 

Consequently, to first infinitesimal order terms the solution /W = + f^ 1 ' coincides with 
the solution /M = + received in the Enskog- Chapman theory: 



(155) 
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As a result, with the same exactness expressions for heat flux density vector (Q, Chapter 7, 
§4) 

qW = q(D i q W = q« = J l - m (c - u™) 2 (c - u™) /Wrfc 



3m dr 

i 2^rwarw 

3m <9r 
dr 



A ■ I (A) dc 
A ■ 7 (a) dc 



-A-r- (156) 



and pressure tensor 

pW ~ pW = y m (c - U W) (c - U W) f^dc 
I ( n (o) kT (o) + n (D fcT (o) + n ^kT^) U 



.-kT^^- u(°) / B : 7(B) dc 
5 ar / 



i nMjfcrWu - i fcT [1] ^ uW y B : 7 (b) 7c 



n WfcrWU-2^u[ 1 ], (157) 



<9r 



coincide, cf. with the contrary assertion, for example, in 151 ] . [121 ] . [131 ]. In fl 157ft notation is 
used: for arbitrary second-rank tensor w 

{^) af 3 = \( W ^ + Wpa) (158) 

- corresponding symmetric tensor. 

Having solutions and we can choose: one can solve separately gas-dynamic 
equations systems of the first and the second orders fll27l) - fll29p and fH9l) . r = 1, and 
separately find n^°\ vS°\ and or, substituting function /W = /(°) + 

in the system of transfer equations (or summing systems of equations fll27l) - fll29p and f|49l) . 
r = 1), immediately seek solutions xv- l \ T' 1 ' of the system (j!17p of, generally speaking, 
singularly perturbed differential equations. 

Analogous results can be received for multi-component gas within the Enskog approach. 
These questions together with third order gas-dynamic equations system within the Enskog 
approach will be considered in a next article. 
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VIII. TURBULENCE AS MULTI-COMPONENT GAS DYNAMICS 



As is well known, laminar flow becomes the turbulent flow, when some parameter char- 
acterizing the flow, namely, Reynolds number 

R=^>1. (159) 

In f ll59p . p is the density of gas, u and L are some characteristic macroscopic velocity and 
linear size of the flow, [i is the coefficient of viscosity. Having rewritten H159[) as 

R= E £, (160) 

- cf. with the expression for viscosity tensor in (11571) . the Reynolds number can be treated as 
the ratio of the macroscopic momentum flux, proportional to f^°\ to the viscosity-induced 
microscopic momentum flux, proportional to Roughly speaking, viscosity "aligning" 
the gas molecules according to a Maxwellian distribution at the same mean velocity and 
temperature can "process" the microscopic momentum flux alone. However, if the macro- 
scopic flux outperforms the microscopic, the gas flow, necessarily, comes to be stratified to 
components. The flow stratification to components can be also caused by external factors. 
If 

f(0) 

i?~4— ^oo, (161) 
f(i) v ' 

then that turbulent flow must be described by the gas-dynamic equations system, corre- 
sponding to the first approximation order in the approximate method for solution of kinetic 
Boltzmann equation, i.e. without viscosity and heat conductivity. But the gas-dynamic 
equations system of the first approximation order in the Enskog- Chapman theory can not 
describe turbulent flow with the entropy increase, see (11321) . Using the gas-dynamic equa- 
tions system (|118j) - (jl20p resolves this paradox. 

If gas-dynamic equations do not describe turbulent gas flows, then either something has 
been missed during the transition from the exact solution of the kinetic Boltzmann equation 
to its approximate solution (by the Enskog method) and then to the gas-dynamic equations, 
or the kinetic Boltzmann equation does not describe turbulent gas flows and requires re- 
placement. However the last, i.e. necessity of replacement of the kinetic Boltzmann equation 
on another kinetic equation at transition from gas laminar flow to turbulent gas flow, seems 
ill-founded. 
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The gas dynamics of the components with the velocity distribution functions, close to 
the Maxwell functions of different mean velocities and temperatures, should be described by 
equations (II 18ft - (11201) . From this point of view, the observed stochasticity of the turbulent 
flow is similar to the stochasticity of the Brownian motion. They differ in scale: in the 
Brownian motion that particle moves stochastically, whose mass is comparable to the mass 
of separate gas molecules, while in the turbulent flow that body moves stochastically, whose 
mass is comparable to the mass of separate gas components. In flll9p - fll20l) the integral 
terms (proportional to n,, rij) can be huge, it explains unexpected power of turbulent effects. 
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